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Abstract. Coherent oscillations of the inflaton field at the end of inflation can para- 
metrically excite fermions in much the same way that bosons are created in preheat- 
ing. Although Pauli-blocking prohibits the occupation number of created fermions 
from growing exponentially, fermion production occurs in a manner significantly dif- 
ferent from the expectations of simple perturbation theory. Here, I discuss the na- 
ture of fermion production after inflation and possible applications including the effi- 
cient transfer of inflaton energy and the production of super-massive fermions during 
' fermionic preheating. 

in 

Consider a simple model of chaotic inflation with the potential |A0 4 for an in- 
flaton field <fi coupled to a massless spin-| field ip by a Yukawa interaction. At the 
end of inflation, the inflaton field will oscillate coherently about the minimum of its 
Q-i! effective potential with an initial amplitude O ~ O(0. 1M P ). In perturbation the- 

ory, one treats the homogeneous and quasi-classical inflaton field as a condensate 
of scalar inflaton particles, each of which can decay into a pair of ^-particles. Each 
fermion then carries away half of the energy of a typical inflaton particle, giving a 
spectrum narrowly peaked around the comoving momentum k m 0A2y/\cj) o . The 
inflaton energy is transfered to fermions after O(rs) 3> 1 inflaton oscillations. 

To investigate fermion production non-perturbatively, we are interested in the 
equation of motion for the field operator i/j. Following the usual prescription (see [1] 
and references therein for details), one seeks eigenfunctions of the Dirac equation 
in the presence of the classical time-dependent source, (j)(t). As the inflaton field is 
spatially homogeneous, only the temporal part of the eigenmode obeys a non-trivial 
equation of motion. These modes, X^t), obey an oscillator-type equation with a 
complex frequency that varies periodically with time: 



X^+( K 2 + qf-i^qf)X k = 



Here, the comoving momentum k enters the equation in the combination k 2 = jio 



and the character of the solutions is defined by the parameter q = — (h is the 
dimensionless Yukawa coupling). The background oscillations enter in the form 
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FIGURE 1. The comoving occupation num- 
ber rifc of fermions in X(f> 4 -inflation as a func- 
tion of time (in units of inflaton oscillations) 
for q= ^ = 1(T 4 and 100 and k 2 = 0.18 and 
11.9, respectively. The period of the modulation 
is about 88 and 22 accordingly. 
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FIGURE 2. The comoving occupation num- 
ber of fermions in m^cj) 2 -inflation as a function 
of (scaled) comoving momentum after 50 infla- 
ton oscillations for initial resonance parameter 
q Q = 10 3 . Expansion destroys the details of the 
resonance band and leads to a fermi-sphere of 

width q 1 / 4 ~ g yv/ 4 . 



/(r) = cn [t, -j=) having unit amplitude and a period T = 7.416. Note that we are 
working in scaled conformal field and time variables so the effects of expansion do 
not appear in the equations of motion. We can express the comoving occupation 
number of particles in a given state through the solutions of eq. (1) 

n k (r) = (Qk ~ Q f f) [ \X' k \ 2 + nl\X k f - 2Q k Im(X k X' k *) ], (2) 

where Vt\ = K, 2 +qf 2 . The energy density of created fermions is = ^3 / d 3 k il k n k . 

It turns out that, for all k 2 and q, the solutions of eq. (1) are periodic in time. This 
is shown by the numerical solutions for the comoving occupation number in Fig. (1). 
Furthermore, it is easy to show that the comoving occupation number defined by 
eq. (2) will obey n k < 1 in accordance with the Pauli principle. We see from Fig. (1) 
that, while the occupation number exhibits some high frequency oscillations (period 
< -|), the most interesting behavior occurs over longer periods. If we average the 
occupation number over an inflaton period, n k (r) = ^ J^ T+T ^ dTn k (r), the average 
occupation number is found to obey the simple equation: n k {r) = F k sin 2 v k r. For 
a given theory, i.e. for a given value of the resonance parameter q, F k < 1 is a 
momentum dependent amplitude and v k is a momentum dependent frequency. In 
Fig. (3) the amplitude F k as a function of k 2 is plotted for several values of q. We 
see that the perturbative expectation is only met for q < 10~ 4 . 

In fact, for large q, the fermions are excited up to k 2 ~ y/q. This is the same result 
as for the broad bosonic resonance and can be understood in the same manner. 
The comoving occupation number, eq. (2), is an adiabatic invariant of the mode 
equation (1). The condition that the modes X k evolve non-adiabatically is Q k > Vl? k 




FIGURE 3. The envelope functions Fk show- 
ing the bands offermion resonance excitation in 
\<j> 4 -inflation for q = ^ = 10~ 4 , 10~ 2 , and 1.0 
(the narrowest to broadest band, respectively). 
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FIGURE 4. The log of the period of modu- 
lation -^j, (in units of inflaton oscillations) as 
a function of k 2 for q = ^- = 1CP 4 , 10~ 2 , and 
1.0 for the light, dotted, and heavy curves re- 
spectively. 



which leads to the condition k 2 < ^J~q for non-adiabatic evolution, and thus, particle 
creation. This occurs much more rapidly than one would expect from perturbation 
theory. In Fig. (4) the period of mode oscillations is plotted in units of inflaton 
oscillations. For all q, the bands saturate after only 10 — 100 inflaton oscillations. 

Turning to the energetics for the most interesting case of the broad resonance 
excitation, q ^> 1, we find ~ O.l/i 2 g 1//4 e0, where the inflaton energy is = j\4>l- 
In chaotic |A0 4 -inflation, A ~ 10~ 13 . If the resonance parameter q is large but 
the coupling parameter is small, h < 0.1, only a small fraction of the inflaton's 
energy will be converted into fermions. Although explosive decay will not occur 
for this model, more general theories such as hybrid models can have large enough 
resonance parameters to allow efficient decay to fermions. 

If the fermion field has a small bare mass term or if we consider inflation with 
a |m|0 2 potential, the conformal invariance of the theory will be broken. In this 
case, the occupation number of fermions for § > 1 no longer evolves periodically 
but becomes stochastic, rapidly fluctuating between = 1 and = 0. This 
destroys the well defined resonance bands depicted in Fig. (3) and can lead to the 
production of super-massive fermions of mass < h<p . These fermions can come 
to dominate the energy density of the universe or survive as massive relics. 

As a particular example of the changes brought by expansion, consider m^cj) 2 - 
inflation with a Yukawa coupling to a still massless fermion. In the broad resonance 
case, one gets a sphere in comoving momentum space with average occupation 
nk = \ that expands with time. A snapshot of this sphere is shown in Fig. (2). 
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